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The problems of using the spectral index of radio galaxies in various tests, in particular, in 
selecting distant radio sources are considered. The history of the question of choosing a criterion 
of searching for distant radio galaxies based on the spectral index is presented. For a new catalog 
of 2442 radio galaxies constructed from NED, SDSS, and CATS data, an analytical form of the 
spectral index. redshift relation has been determined for the first time. The spectral index. angular 
size and spectral index. flux density diagrams have also been constructed. Peculiarities of the 
distribution of sources on these diagrams are discussed. 
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1. INTRODUCTION 

Radio galaxies are among the most powerful observed 
space objects, making it possible to use such radio 
sources to investigate the properties of the Universe 
at various cosmological epochs. Therefore, construct- 
ing samples of radio galaxies from various red- shift 
ranges is one of the most important tasks in observa- 
tional cosmology. 

The sharp fall in the number of observed radio 
sources with increasing redshift forces us to seek for 
ways of rapidly selecting distant objects. One of the 
ways is selection by the spectral index a < —1.0 
(S ~ v a , where S is the flux density and v is the 
frequency). It is based on the fact that the farther 
the object, the more probable that it will have a 
steep spectrum. This is one of the first discovered 
and strongest criteria in searching for distant galaxies 
based on radio-astronomical data. It was established 
independently in several works devoted to the iden- 
tification of radio sources and the analysis of radio- 
spectra statistics. Among the first papers, we will note 
the paper by Whitfield (1957), who pointed out a cor- 
relation of the spectral index for a radio source with 
its distance, and the papers by Dagkesamanskii (1969, 
1970), who found that there were no distant objects 
with spectral indices a > —0.7 for the 3C sources. 
Subsequently, Tielens et al. (1979) determined that 
the fraction of optically identified sources with ul- 
trasteep spectra (spectral indices af^ < —1) de- 
creased with decreasing a. Following this work, Blu- 
menthal and Miley (1979) showed the spectral index 
for 3C and 4C radio sources of various populations 
to depend on the properties of the objects: their ap- 



parent magnitude, redshift, radio luminosity, and an- 
gular size. The detected correlation suggested that 
the steep-spectrum sources were, on average, further 
away and more luminous than the sources with less 
steep spectra of the same populations (radio galax- 
ies and quasars). Laing and Peacock (1980) investi- 
gated the relation between radio spectrum and ra- 
dio luminosity for samples of extragalactic sources at 
178 and 2700 MHz. The spectra were measured for 
the extended regions of the radio sources, which were 
classified by morphological types. At low frequencies, 
the degree of spectral curvature was found to corre- 
late with the luminosity for sources with hot spots. 
At high frequencies, the correlation between spectral 
index and luminosity was confirmed. Laing and Pea- 
cock (1980) also confirmed the existence of a spec- 
tral index. redshift relation at 178 and 2700 MHz for 
FR II sources (Fanaroff and Riley 1974). At present, 
many groups and authors use such a selection of can- 
didates for distant radio galaxies (see, e.g., Cham- 
bers et al. 1988; Wieringa and Katgert 1991; Soboleva 
and Temirova 1991; Rottgering et al. 1997; de Breuck 
et al. 2000; Pedani 2003; Verkhodanov et al. 2003; 
Copal-Krishna et al. 2005; Klamer et al. 2006; Bor- 
nancini et al. 2006; Kopylov et al. 2006). Although 
the criterion based on the spectral index is very effi- 
cient, its explanation is not yet completely clear (De 
Young 2002). Three main, widely used ideas that ex- 
plain this dependence can be highlighted: 

- if the spectral steepness for many radio galaxies 
increases with frequency, then the spectra of distant 
objects must be steeper than those of near ones be- 
cause of the factor (1+z); 
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- the losses due to the Compton scattering of cos- 
mic microwave background (CMB) photons by rela- 
tivistic electrons grow, since the CMB radiation den- 
sity increases as (1 + z) 4 ; these growing losses lead to 
an aging of the population of electrons, causing a cut- 
off or an increase in the steepness of the high-energy 
part of the spectrum where the losses are greatest; 

- the selection effect: only bright sources are seen 
at high redshifts, and it can then be said that pre- 
cisely these sources exhibit the greatest "depletion" 
in the high-energy part of the spectrum. 

Klamer et al. (2006) give yet another possible ex- 
planation for the a — z correlation. They noted that 
the steep-spectrum sources are rather rare among the 
nearest FR I radio galaxies, while those with such 
spectra are located in regions with a high baryon 
density. It can then be assumed that if there is evo- 
lution of the environment of powerful radio galaxies 
reflected in the richness of the surrounding cluster, 
then, on average, the radio galaxies are more likely 
located in regions with a higher ambient density than 
in less dense regions. Hence follows the observed red- 
shift dependence of the spectral steepness. This can 
play its role when the gas density as one goes into 
the past increases as (1 -I- z) 3 , while the injection of 
electrons with steeper spectra occurs naturally and 
is described as a function of the redshift in the first- 
order Fermi acceleration processes attributable to the 
decreasing expansion velocity of the hot spots in the 
denser and hotter intergalactic medium (Athreya and 
Kapahi 1998). 

Other authors (Kharb et al. 2008), who studied 
a sample of 13 "powerful classical double" FR II ra- 
dio galaxies, also analyzed this a — z correlation and 
pointed out the existence of a redshift dependence of 
the spectral index for both the hot spots and the core. 
From the presence of a correlation between the spec- 
tral indices for the hot spots and a core with a flatter 
spectrum, Kharb et al. (2007) drew the conclusion in 
favor of choosing the explanation of the effect by the 
electron aging model. 

Although as yet there is no general, common con- 
sensus on clarifying the nature of the a— z correlation, 
the very fact that this empirical dependence exists is 
important for studies. Almost all of the distant radio 
galaxies found has passed through the stage of such 
a selection. An example of an object selected in this 
way is the radio galaxy that is the red- shift record- 
holder with z=5.19 and spectral index a = —1.63 
(van Breugel et al. 1999). Another example is the 
object RC 0311+0507 investigated in the "Big Trio" 
Program (Kopylov et al. 2006). Its redshift 

is z=4.514, it has the record radio luminosity for 
radio galaxies with z > 4, and its spectral index is 
a = -1.33. 



We prepared a sample of distant (z > 0.3)radio 
galaxies (Khabibullina and Verkhodanov 2009a, 
2009b, 2009c) using the NEE0 CATS^ (Verkhodanov 
et al. 2005), and SDSS^ (Schneider et al. 2007) 
databases for the subsequent application in various 
statistical and cosmological tests (Verkhodanov and 
Parijskij 2003, 2009; Miley and De Breuck 2008), in 
which a large number of objects of the same nature 
is required to carry out a study. To compile the pri- 
mary list, we used the NED database, from which we 
selected objects with the following parameters: the 
redshift (z > 0.3) and morphological properties. radio 
galaxies. The initial list contained 3364 object. This 
sample of objects is contaminated by objects with 
incomplete information or by objects with different 
properties. Therefore, the next stage was to clean the 
initial sample of superfluous sources. For this purpose, 
we selected the objects that were removed from the 
primary list (Khabibullina and Verkhodanov 2009a): 
(1) with photometrically determined redshifts and (2) 
with quasar properties based on available published 
data. The final catalog contains 2442 sources with 
spectroscopic redshifts, photometric magnitudes, ra- 
dio flux densities, sizes of the radio sources, and ra- 
dio spectral indices calculated from the results of the 
cross-identification of the list of selected radio galax- 
ies with the radio catalogs stored in CATS in the fre- 
quency range from 30 GHz to 325 MHz. By default, 
the flux densities are given on the Baars scale. The 
data collected in the catalog can be used for cosmo- 
logical luminosity- redshift, angular size- redshift, and 
age-redshift tests. 

In addition, the statistical diagrams for the pa- 
rameters of radio galaxies and their evolutionary 
properties can be investigated using the cataloged 
data. 

The goal of this paper is to construct and analyze 
the spectral index-redshift, spectral index-angular 
size, and spectral index-flux density diagrams for a 
pure sample of radio galaxies, which can be used to 
estimate the redshifts of radio galaxies and to calcu- 
late their luminosity function. 

2. CALCULATING THE SPECTRAL 
INDEX 

To calculate the spectral indices, we performed cross- 
identification in the CATS database with a 200" x 
200" identification window. To remove chance field 
radio objects in the specified box, we used a technique 



1 http://nedwww.ipac.caltech.edu 

2 http://cats.sao.ru 

3 http://www.sdss.org 
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of data analysis similar to that described 
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Figure 1: Histogram of the spectral-index distribution 
at 1400 GHz. 

in Verkhodanov et al. (2000, 2009). The essence 
of the method is to use a joint analysis of the data 
in coordinate and spectral spaces to separate out the 
probable identifications of specific radio sources at 
various radio frequencies. The spg program (Verkho- 
danov 1997a) of the RATAN-600 continuum data 
processing system is used for these purposes. In de- 
scribing the spectra S(u) of sources with available 
measurements at several frequencies, for the subse- 
quent calculation of the spectral indices we used an 
automatic parametrization of S(v) by the formula 
logS{v) = A + Bx + Cf(x), where S is the flux den- 
sity in Jy, x is the logarithm of frequency v in MHz, 
and f(x) is one of the following functions: exp(-x), 
exp(x),or a; 2 . For 91our sample, the spectra are linear. 
The accuracy of determining the coefficients of the 
least-squares fits varies from spectrum to spectrum, 
but it does not exceed 10a given frequency was deter- 
mined as the slope of the tangent to the spectrum at 
the specified point. 

Figure 1 presents a histogram of the spectral- in- 
dex distribution at 1400 MHz. Since the spectra of the 
sample were fitted mainly by linear functions, the dis- 
tribution is virtually independent of the wavelength. 
The median of the spectral-index distribution for our 
entire sample is -0.63. Figure 2 presents the redshift 
distribution of objects. 

3. STATISTICAL DIAGRAMS 

Using the list of selected radio galaxies (Khabibullina 
and Verkhodanov 2009a) , we constructed the spectral 
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Figure 2: Histogram of the redshift distribution for 
the catalog. 

index-redshift (a— z), spectral index-angular size (a — 
6), and spectral index- flux density (a — S) diagrams. 

The spectral index-redshift diagram is shown in 
Fig. 3. It has a clear tendency for the spectral index 
to decrease with increasing z and satisfies a linear 
regression a = a + 6z,where a = —0.7 ± 0.02 is the 
regression constant and b = —0.15 ± 0.01 the slope. 
Some of the objects have peaks in their spectra at a 
frequency near 1 GHz and lower, giving a 
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Figure 3: Spectral index-redshift diagram. The spec- 
tral indices were calculated at 1400 MHz. The trian- 
gles and circles indicate the SDSS objects and the ob- 
jects from all the remaining catalogs, respectively. The 
regression line was drawn using the median spectral 
indices within bins with a step Az = 0.5. 

positive spectral index at low frequencies. In the 
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Figure 4: Pseudo-three-dimensional angular size- 
redshift- spectral index diagram for radio galaxies. 



case of negative a, when the source has a steep spec- 
trum, a straight line was used to fitit. The diagram 
was constructed for the spectral indices at 1400 MHz. 
The presented boundaries of the parameter estimates 
correspond to the statistical ones determined by the 
least-squares method from the median values. The ac- 
curacy of determining the regression parameters can 
be discussed in terms of sample completeness. How- 
ever, as the above estimates show, the extrapolation 
prediction gives a difference in the right boundary 
value of less than 13from the first half of the bins. 
The selection of distant objects by the steepness of 
their spectra has been widely used in a number of 
papers (Blumenthal and Miley 1979; Parijskij et al. 
1996, 1999; de Breuck et al. 2000), but the analyti- 
cal form of the regression for such a large sample of 
radio galaxies has been determined for the first time. 
The regression was calculated from the median spec- 
tral indices in zintervals with a step (bin) Az = 0.5. 
The choice of the bin size Az — 0.5 was dictated by 
two factors: it should not be (1) large, Az < 1.0, in 
order to take into account the differential peculiari- 
ties in the behavior of the population and (2) small, 
Az > 0.1, in order that there be a sufficient number 
of objects in the specified bins. We chose Az = 0.5, on 
the onehand, as the mean size and, on the other hand, 
as a bin that allowed a comparatively large number of 
objects to be obtained at z > 3. The chosen bin size 
makes the regression parameters very stable against 
an increase in the number of objects during a further 
possible expansion of the sample. 

In principle, the existence of a slope can be ex- 
plained by two selection effects, although it docs not 



rule out physical causes either. 
— (1) The lists of distant radio galaxies were taken 
from the catalogs that were compiled by taking into 
2 s account the selection by a steep spectrum. This effect 
2 can be avoided only based on complete samples con- 
1 taining radio galaxies with measured redshifts. There 
are very few such samples. (2) When we choose the 
most powerful radio sources at high z, hot spots with 
steep spectra contribute to the radio emission. The 
farther the radio source, the more probable that it 
will be more powerful than the surrounding ones and 
will have a steep radio spectrum. We are planning 
to verify this fact in our subsequent paper by study- 
ing subsamples of galaxies with different luminosities. 
The analytical form of the regression can also be used 
to preselect objects at given z, to estimate the dis- 
tances to radio galaxies in the first approximation, 
and to analyze and model the distributions of the 
number of sources (Condon 1984; Gorshkov, 1991), 
which are used to study the luminosity function. The 
estimate of the regression parameters is stable, be- 
cause it was made using the median values of large 
subsamples of radio galaxies. It should be noted that 
when the spectral indices of radio sources from sam- 
ples including both radio galaxies and quasars are 
analyzed, the regions of .with low zare filled (Vollmer 
et al. 2005). This shifts greatly the dependence a(z) 
toward larger spectral indices. In addition, to check 
the stability of the slope, we limited the data used 
by redshift intervals with z < 2.5 (where more than 
15 objects fell within the bins being analyzed) and 
reconstructed the observed regression for them. As a 
result, the value of .predicted by the regression with 
incomplete data for z=5.2 was -1.34, in contrast to 
-1.54 obtained for the complete set of bins, i.e., the 
relative error in the extrapolation to high zis within 
13consider satisfactory. 

For the angular size. redshift. spectral index dia- 
gram (Fig. 4), we made the spectral index. angular 
size (Fig. 5) and spectral index. flux density (Fig. 6) 
projections; wc presented the former in two forms, 
highlighting in one case the SDSS objects by symbols 
and in the other case the objects from different red- 
shift ranges by different symbols: 0.3 < z < 0.7 (cir- 
cles), 0.7 < z < 1.5 (squares), and 1.5 < z (asterisks). 
The angular sizes were taken from the NVSS cata- 
log (Condon et al. 1998). The last diagram demon- 
strates the existing relation between the red- shift of 
a radio galaxy and the radio spectral index. In addi- 
tion, the expected relation where larger sources have 
steeper spectra due to synchrotron radiation in ex- 
tended components is observed in Fig. 5. Note that for 
the distribution of radio galaxies observed in SDSS, 
there is a shift toward smaller spectral indices for the 
case where the distributions of sources in size for the 
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Figure 5: (a) Spectral index- angular size diagram. 
The triangles and circles indicate the SDSS objects 
and the objects from all the remaining catalogs, re- 
spectively, (b) Spectral index-angular size diagram. 
Differen symbols indicate different redshift ranges: 
0.3 < z < 0.7 (circles), 0.7 < z < 1.5 (squares), 
and 1.5 < z (asterisks). The spectral index was calcu- 
lated at 14-00 MHz. The measured sizes of the radio 
sources were taken from the NVSS catalog. 



SDSS and non-SDSS subsamples do not differ in form. 
The difference in the data shown in Fig. 5 can be ex- 
plained by the selection effect attributable to the ex- 
istence of a spectral index. redshift relation. The latter 
is demonstrated in Fig. 5b, where the region of high 
zis marked by light gray. 

On the spectral index-flux density diagram (Fig. 
6), we marked galaxies with redshifts in different 
ranges by different symbols: 0.3 < z < 0.7 (circles), 
0.7 < z < 1.5(squares), and 1.5 < z asterisks). The 
flux densities were taken from the NVSS catalog (1400 
MHz). As we see from the figure, objects from all 
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Figure 6: Spectral index-flux density diagram for ra- 
dio galaxies. The flux density was taken from the 
NVSS catalog at 1400 MHz. Different symbols mark 
objects from different redshift ranges: 0.3 < z < 0.7 
(circles), 0.7 < z < 1.5 (squares), and 1.5 < z (aster- 
isks). 



redshift ranges are encountered in the entire range 
of flux densities; most of the objects with large spec- 
tral indices (a > —0.3) are the objects with the low- 
est redshifts in our sample. In addition, a correlation 
<*('Si40o) between the spectral index and flux density 
is observed for all three zranges. We drew linear re- 
gressions in the form a — d + klogSuo$ using the 
median values of a calculated in logarithmic intervals 
of flux densities. The regressions for different ranges 
are indicated in Fig. 6 by the solid line (d = —0.60; 
It = -1.8 x 10~ 4 ) for 0.3 < z < 0.7, the dashed line 
(d = -1.03; k = -3.7 x 10~ 5 ) for 0.7 < z < 1.5, 
and the dotted line (d = -1.04; k = -9.2 x 10~ 5 ) for 
1.5 < z. They confirm the apparent effect of a change 
in the slopes of the correlations for different redshifts. 



4. CONCLUSIONS 

We presented the results of our investigation of the 
a — z, a — 0, and a — S diagrams for a sample of radio 
galaxies containing 2442 objects (Khabibullina and 
Verkhodanov 2009a, 2009b, 2009c) and constructed 
from the spectroscopic data collected in NED, SDSS, 
and archives of the Big Trio program (Parijskij et al. 
1996, 1999). For the objects of our list, we found the 
existence of a regression ct(z) and established its an- 
alytical form: a(z) — —0.73 — 0.15z. The estimate of 
the regression parameters is stable, since is was made 
using the median values of large subsamples of radio 
galaxies. This dependence can also be used to pres- 
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elect objects at given z,to estimate the distances to 
radio galaxies in the first approximation, and to study 
the luminosity function. Nevertheless, when using the 
derived relation ot(z), we should take into account the 
fact that the detected regression was obtained from 
an incomplete sample. 

For the correlation a(S'i4oo)j the regression pa- 
rameters were found to evolve with increasing z. Since 
correlations (Disney et al. 2008) suggesting the exis- 
tence of a multiparameter fundamental plane of ra- 
dio galaxies are found for a number of their physical 
parameters (total mass, baryon fraction, luminosity, 
etc.), it may be concluded that the evolution of the 
dependence a(5i4oo) is also indicative of the possi- 
ble evolution of the fundamental plane. In turn, the 
existence of a correlation between the variations in 
parameters of the fundamental plane and spectral in- 
dices argues for the explanation of the dependence 
a(z) by the electron aging model, which complements 
the conclusions by Kharb et al. (2007). 
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